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ABSTRACT 

We present results from a study of a nuclear emission of a nearby radio galaxy, 4C+29.30, over 
a broad 0.5-200 keV X-ray band. This study used new XMM-Newton (~ 17ksec) and Chandra 
(~ 300ksec) data, and archival Swift/BAT data from the 58- month catalog. The hard (> 2keV) X- 
ray spectrum of 4C+29.30 can be decomposed into an intrinsic hard power-law (r ~ 1.56) modified by 
a cold absorber with an intrinsic column density Ah. z ~ 5x 10 23 cm -2 , and its reflection (|f2/27r| ~ 0.3) 
from a neutral matter including a narrow iron Ka emission line at the rest frame energy ~ 6.4 keV. 
The reflected component is less absorbed than the intrinsic one with an upper limit on the absorbing 



column of Af \ < 2.5 x 10 



The X-ray spectrum varied between the XMM-Newton and 



Chandra observations. We show that a scenario invoking variations of the normalization of the power- 
law is favored over a model with variable intrinsic column density. X-rays in the 0.5-2 keV band are 
dominated by diffuse emission modeled with a thermal bremsstrahlung component with temperature 
~ 0.7 keV, and contain only a marginal contribution from the scattered power-law component. We 
hypothesize that 4C+29.30 belongs to a class of 'hidden' AGN containing a geometrically thick torus. 
However, unlike the majority of them, 4C+29.30 is radio-loud. Correlations between the scattering 
fraction and Eddington luminosity ratio, and the one between black hole mass and stellar velocity 
dispersion, imply that 4C+29.30 hosts a black hole with ~ 10 s M mass. 

Subject headings: galaxies: active — galaxies: individual (4C+29.30) — X-rays: galaxies — accretion, 
accretion disks 
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1. INTRODUCTION 

4C+29.30 is a low-redshift (z = 0.0647) radio 
source with a moderate radio luminosity ~ 10 42 ergs _1 
(Siemiginowska et al. 2012) hosted by an elliptical 
galaxy. It was first studied in the radio and optical 
bands by van Breugel et al. (1986). Observed complex 
radio morphology was resolved into jets, knots, lobes, 
and a diffuse tail. A number of optical emission lines 
were measured, and their relative intensities indicated 
photoionization by a non-stellar continuum attributed 
to the central active galactic nucleus (AGN). The ra- 
dio emission and the optical line emitting gas were found 
to be strongly linked suggesting interaction between the 
radio source and its environment. In addition, shells and 
dust lanes observed in the optical band (van Breugel 
et al. 1986) provide evidence of a past merger with a 
gas-rich disk galaxy which possibly triggered the rejuve- 
nated AGN activity. Subsequently, follow-up radio stud- 
ies revealed radio structures with diverse spectral ages 
suggesting an intermittent nature of the radio source in 
4C+29.30. Jamrozy et al. (2007) presented evidence for 
a large scale extended relic radio emission most probably 
due to an earlier cycle of activity of the source > 200 Myr 
ago. The age of a small-scale radio structure embedded 
in the extended relic radio emission was estimated at 
< 100 Myr, with the inner double knots of a spectral age 
of < 33 Myr (Jamrozy et al. 2007). Liuzzo et al. (2009) 
resolved the very central region of the source into two 
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nuclear knots with spectral ages of ~ 15 yr and ~ 70 yr. 
4C+29.30 was cataloged as an infrared point-source by 
the IRAS, WISE and 2MASS surveys (NASA/IPAC In- 
frared Science Archive). 

The first detection of the source in the X-ray band 
occurred during a snap-shot 8 ksec Chandra observation 
(Gambill et al. 2003). The short exposure time and 
few detected counts prevented detailed spectral analy- 
sis. Nevertheless, this pilot observation hinted towards 
a complex X-ray morphology composed of a nucleus, 
hotspots, a putative jet, and diffuse emission. The 
source was not a target of X-ray pointings prior to the 
first Chandra observation, without detection in, e.g., the 
ROSAT All Sky Survey, probably due to significant in- 
trinsic absorption (Gambill et al. 2003, Siemiginowska et 
al. 2012). In 2008, we observed 4C+29.30 with XMM- 
Newton for ~ 45 ksec. This observation was followed 
by our deep 300 ksec Chandra exposure (Siemiginowska 
et al. 2012, hereafter Paper I) which clearly revealed 
diverse X-ray emitting structures, many of which cor- 
respond to those in radio or optical bands, and others 
are intrinsic to the X-ray band. The soft 0.5-2 keV X- 
ray image showed northern and southern lobes, hotspots, 
southern jet, and thermal diffuse emission in the center. 
Emission in the hard X-ray band (above 2 keV) was in- 
stead dominated by the active nucleus of 4C+29.30. In 
Paper I, our discussion focused on the spectral proper- 
ties of the soft X-ray components, the structure of the 
outflow, and interactions of the X-ray jet with the in- 
terstellar medium. We emphasized similarities between 
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4C+29.30 and famous radio galaxies showing jet-ambient 
medium interactions, i.e. M87 (e.g., Million et al. 2010), 
NGC 1275 (e.g., Fabian et al. 2011), and Cen A (e.g., 
Morganti et al. 2010, and references therein). 

In the present paper, we focus on detailed modeling 
of the X-ray core emission to uncover the nature of the 
AGN powering 4C+29.30. We analyze new Chandra data 
and archival Swift/BAT data from the 58-months cata- 
log (Baumgartner et al. 2010), and thus we study the 
core emission over a broad 0.5-200 keV X-ray band. We 
apply our best fit model to our new XMM-Newton obser- 
vations of 4C+29.30 taken ~ 2 years before our Chandra 
observation. We search for signatures of the X-ray spec- 
tral variability between the Chandra and XMM-Newton 
observations. The paper is organized as follows. In Sec- 
tion 2 we describe our data and provide details on the 
preparation of the spectra for modeling. In Section 3 we 
test various spectral scenarios for AGN X-ray activity. 
We compare the best fitting Chandra and XMM-Newton 
models in order to study the origins of possible spectral 
variability. In Section 4 we present the discussion of our 
results indicating that 4C+29.30 is a strongly absorbed 
intrinsically X-ray variable AGN with Seyfert 2 type of 
activity and a peculiar obscuring geometry, probably due 
to a geometrically thick torus similar to that reported for 
several Swift/BAT selected AGN (e.g., Ueda et al. 2007; 
Winter et al. 2009a, b). Finally, in Section 5 we summa- 
rize and conclude our study. Throughout the paper we 
adopt Q A = 0.73, n M = 0.27 and H = 70kms _1 Mpc -1 
for the flux-luminosity transformations. 

2. X-RAY OBSERVATIONS 

2.1. XMM-Newton 

Our XMM-Newton observation of 4C+29.30 was per- 
formed on Apr 11, 2008. It was split in two pointings 
lasting - 45.7 ksec (ID: 0504120101) and - 9.7ksec (ID: 
0504120201). No science data were produced during the 
second pointing. The data reduction of the 45.7 ksec long 
observation was performed following the standard proce- 
dure!]] with the XMM-SAS v. 11.0.0 package using the up- 
dated calibration files (Oct 2011). We constructed light 
curves at energies above 10 keV to identify periods of high 
background flaring. After cleaning the data, we obtained 
the following net exposures and 0.5-10 keV observed 
count rates: ~ 17ksec and (6.3 ± 0.2) x 10 _2 cnts _1 for 
the EPIC/PN camera, - 23 ksec and (1.75±0.09) x 10~ 2 
cnt s- 1 for MOS1, and - 23 ksec and (2.03 + 0.09) x 10" 2 
cnt s _1 for MOS2. In this paper we will limit the discus- 
sion to the EPIC/PN data which provide better statistic 
than the MOS data. However, the data taken by the 
EPIC/MOS1 and MOS2 have been checked for consis- 
tency. 

The source 0.5-10 keV PN spectrum was extracted 
from a circular region (centered at R.A. 08:40:02.276, 
Decl. +29:49:02.47; J2000.0) of 25" radius (Figure [TJ, 
and the background from a circular region located on the 
same CCD (50" radius; not indicated in Figured]). The 
response matrices were created using the SAS commands 
rmf gen and arf gen. The source data are not piled-up. 
We used specgroup to set spectral groups such that each 
group/bin contains a minimum signal-to- noise ratio of 3. 

1 http:/ /xmm.esac.esa.int/sas/current/documentation/threads/ 



2.2. Chandra 

Deep Chandra ACIS-S imaging observations of 
4C+29.30 were performed in Feb 2010, i.e. ~ 2 years 
after our XMM-Newton observation. The approved 
300 ksec observation was split into four separate point- 
ings (OBSIDs 11688, 11689, 12106, 12119) that sum to a 
total exposure of 286.4 ksec (see Paper I for details of the 
observations and data reduction). In this paper we use 
the four combined^ 0.5-7 keV X-ray spectra (5328 ± 73 
total counts) of the core region (Figure QJ defined as a 
circle with radius of 1.25" centered at the position of 
the nucleus determined based on the hard band Chan- 
dra X-ray image (R.A. 08:40:02.345, Decl. +29:49:02.61, 
J2000.0; Paper I), together with the appropriate calibra- 
tion response files. The background (25 ± 5 counts) was 
extracted from an annulus centered on the core position 
with an inner radius of 1.5" and an outer radius of 10". 
We use CI AO version 4.3 tools for the data analysis. 

2.3. Swift/BAT 

We extend the Chandra and XMM-Newton energy 
bandpasses by considering the publicly available 14- 
195 keV (8-channel) hard X-ray Swift/BAT spectrum 
from the 58-months catalog (between Nov 2004 and Sep 
2008; Baumgartner et al. 2010) together with the corre- 
sponding response matrix diagonal_8 . rsp. We do not 
detect any long term variability signal in the Swift/BAT 
lightcurve (normalized excess variance cr 2 xs < 0; e.g., 
Vaughan et al. 2003). Either way, integrating the data 
over the period of ~ 5 years would smear any hard X- 
ray spectral variations. Thus, the Swift/BAT data are 
representative of the average hard X-ray emission from 
4C+29.30. 

3. MODELING 

In this section we use our deep Chandra observation 
to establish a spectral model for the 4C+29.30 core con- 
tinuum. We use the Swift/BAT data to check the con- 
sistency of the model with X-ray data outside the Chan- 
dra energy bandpass. Finally, we apply the best fitting 
Chandra model to the XMM-Newton observation taken 
~ 2 years prior to the Chandra pointing, and we discuss 
the X-ray spectral variability of the source. We perform 
the fits in Sherpa (Freeman et al. 2001) using CSTAT 
statistic. 

The Swift/BAT data were averaged over 58 months 
of observations, which potentially might introduce a 
discrepancy between the Chandra (or XMM-Newton) 
and Swift/BAT data due to the X-ray flux variability. 
Swift/BAT composite spectra contain fractional counts 
with Gaussian statistic and are background subtracted, 
which prevents the CSTAT fitting. For this reason, 
we find best fit of each model to the Chandra (or XMM- 
Newton) data only. Then we apply the best fit model 
to the Swift/BAT data using % 2 statistic with variance 
supplied with the data chi2xspecvar, allowing for only 
one free parameter, the normalization constant, Cbat- 

2 Because the data were collected on the same part of the ACIS- 
S detector and spanned a relatively short time interval, we used 
CIAO 4.3 script combine_spectra to combine individual spectra 
and response files into one spectrum for this analysis. 

3 http: / /heasarc. nasa.gov/docs/swift/analysis/threads/batspectrumthread.ht: 



Nuclear X-ray properties of the peculiar radio-loud hidden AGN 4C+29.30 



3 



Table Q] contains the description of all models used in 
the paper in the XSPEC (Arnaud 1996) and Sherpa termi- 
nology. In Table [2] we list the best fit model parameters 
and the resulting fit statistics. All models considered in 
this section include Galactic absorption (phabs) with a 
column density fixed at Ah = 3.98 x 10 20 cm" 2 (Dickey 
& Lockman 1990). Errors on the spectral parameters 
correspond to the 68% confidence level (la) for one sig- 
nificant parameter. 

A region of the diffuse emission surrounding the nu- 
cleus can be clearly identified in the soft 0.5-2 keV Chan- 
dra image ('Center' in Figure 2 in Paper I). In Paper 
I, we showed that a thermal bremsstrahlung model pro- 
vided a better fit to the spectra extracted from this region 
than a power-law model. Thus, all our models include 
a thermal bremsstrahlung component with the tempera- 
ture and normalization allowed to vary. 

We modeled the Chandra and XMM-Newton back- 
ground spectra by fitting appropriate background regions 
(Figure [1]) with a model consisting of a Galactic ab- 
sorption, thermal bremsstrahlung and power-law. The 
background model parameters are listed in Tabled The 
parameters differ because the XMM-Newton background 
was collected from a region at the distance > 25" from 
the nucleus. 

3.1. Chandra 

We considered several models for the nuclear X-ray 
emission of 4C+29.30 observed with Chandra. 

Model 1. First, we tried a simple parametrization us- 
ing a cut-off power-law function (zhighect*zpowerlw in 
XSPEC terminology; in the following we will refer to this 
component as CPL or 'a power-law'). The zhighect 
model is defined as exp{E\ / E{ a \& — E x (1 + z)/£foid]- 
We assumed E\ — 0. The E-folding energy, £f id could 
not be constrained by our data and we fixed it at 300 keV. 
We added an intrinsic absorber with a variable column 
density (zphabs), motivated by the presence of a dust 
line crossing the nuclear region, as revealed by the HST 
image (Figure^. We were not able to obtain a good 
fit with this model. The plot of the model-to-data ra- 
tio revealed residuals in the 5-7 keV band suggesting the 
presence of an iron Ka emission line (Figure |3^,) . These 
residuals forced the photon index of the intrinsically ab- 
sorbed (Ah, z ~ 4 x 10 23 cm -2 ) continuum to converge 
to an unusually low value, T = 0.84+°; 27. 

Model 2. In Model 2 we added a Gaussian line 
(zgauss) which reduced the 5-7 keV residuals signif- 
icantly (Figure [3}d)- The resulting line was narrow, 
a = (0.07 ± 0.03) keV, and centered at the rest frame 
energy of Ek u = (6.36 ± 0.02) keV. These parameters 
suggest that the line originates in a neutral material far 
from the black hole, e.g., at the outskirts of an accre- 
tion disk, or a dusty torus. The earlier short 8ksec 
Chandra observation contained too few counts to detect 
the line or put an upper limit on its EW (Gambill et 
al. 2003). Thus, this is the first detection of the iron 
Ka line in 4C+29.30. The observed equivalent width 
of the line (i.e. relative to the intrinsically absorbed 
continuum) was EW = 192+ig eV, while the intrinsic 
one (i.e. relative to the unabsorbed continuum) was 
EW; = 9li 2 geV. The photon index in Model 2 was still 
unusually low, T = 0.97tom, 

and the intrinsic absorp- 



tion remained strong with Ah, z ~ 5 x 10 23 cm~ 2 . The 
temperature of the diffuse emission, kT = 5.2+ 3 'gkcV, 
was in agreement with that measured in Model 1 (kT — 
6.4+2'ikeV), but considerably higher than the tempera- 
ture of the background thermal bremsstrahlung compo- 
nent, kT = (0.63 ± 0.04) keV. 

Model 3. The detection of the iron Ka line suggested 
that the total spectrum may contain a certain contri- 
bution of the associated reflected emission. Hence, we 
considered the reflection model with a slab geometry, 
pexriv of Magdziarz & Zdziarski (1995), i.e. a reflec- 
tion from an accretion disk. We first assumed that the 
reflected component is not intrinsically absorbed, con- 
trary to the case of the incident CPL. This corresponds 
to a scenario in which the CPL continuum is observed 
through an obscuring material, e.g., a dusty torus, while 
the reflected emission arrives to the observer unobscured. 
With default settings the pexriv model computes both 
the incident CPL and the reflected component. However, 
in order to apply the intrinsic absorber only to the inci- 
dent CPL, we leave the CPL component in the model, 
and choose a setting in which the pexriv model returns 
only the reflected component. Consequently, the reflec- 
tion amplitude, Q/2tt, was fitted over a negative range, 
which physically corresponds to a reflection component 
with amplitude |f2/27r|. 

We linked the photon index and the normalization of 
the pexriv model to those of CPL. The pexriv reflec- 
tion does not account for the iron line emission and so we 
kept the zgauss component explicitly in the model. The 
viewing angle, i, is not known for our source. However, 
an upper limit can be calculated based on the jet-to- 
counterjet flux ratio, J, at 5 GHz (e.g., Urry & Padovani 
1995). Sambruna et al. (2004) and Liuzzo et al. (2009) 
give J > 29 and J > 50, implying i < 60° (for dis- 
crete knots). We fixed i at 45°, assuming that the jet in 
4C+29.30 is well aligned with the axis of the accretion 
disk. We fixed the abundances of elements at the Solar 
value. 

The resulting data-to-model residuals are presented 
in Figure [3fc. The ionization parameter in the pexriv 
model converged to zero. The iron line remained nar- 
row, and centered at the rest frame energy of E^ a — 
(6.36 ± 0.01) keV. Its observed EW„ = mtgeV (in- 
trinsic EW; = 71± 2 jeV) was lower than in Model 2. 
Thus, the properties of the iron line and the reflected 
emission support the hypothesis that both components 
originate in a neutral matter located far away from the 
black hole. The power-law component with the pho- 
ton index of T = 1.561^1 was absorbed with an in- 
trinsic column density Ah, z ~ 5 x 10 23 cm~ 2 . The 
reflection amplitude yielded an intermediate value of 
|fi/27r| = 0.32±g;?|. The temperature of the soft dif- 
fuse emission, kT = 0.71+0 12 keV, was significantly lower 
than that in Models 1 and 2 due to the contribution of the 
reflected component to the soft X-ray band, and slightly 
higher than (but consistent with) the temperature of the 
background thermal component. As a final model check 
we replaced the thermal bremsstrahlung component with 
an apec model. We were able to establish only an upper 
limit on the value of metal abundances, A < 0.01. It is 
generally known that low abundances might be an arti- 
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fact resulting from fitting the mixed temperature plasma 
where these mixed temperature components cannot be 
spectrally or spatially resolved (see e.g. Kim 2012 in 
context of abundances in hot ISM of elliptical galaxies). 
This problem with regard to the data of 4C+29.30 was 
considered in more details in Siemiginowska et al. (2012). 
Similarly to their approach we decided to keep the ther- 
mal bremsstrahlung component in the model as a good 
approximation. 

Model 4- Next, we allowed for the second intrinsic ab- 
sorber altering the reflected component with a column 
density A^, in addition to the one modifying the in- 
trinsic CPL emission. This model, however, converged 
to Model 3 with N£ R Z = 1-2±^;| x 10 22 cm~ 2 , consistent 
with zero. 

Model 5. Motivated by the studies of other Seyfert 2 
type galaxies (e.g., Evans et al. 2007, 2010; Grandi et al. 
2007) we tested if a partial covering model (representa- 
tive for, e.g., a patchy torus) could provide an alternative 
description of our data. We assumed that only a fraction, 
f pc (allowed to vary between and 1), of the power-law 
emission is passing through an obscuring matter, while 
the remaining fraction (l-/ pc ) arrives to the observer un- 
absorbed. This model also converged to Model 3 with 

/ pc - 1.0000-0.0021. 

Model 6. Next, we checked if a certain fraction, / sc , of 
the intrinsic CPL may contribute to the soft emission be- 
low 1 keV due to scattering. We found that the scattered 
component was not formally required by the data with an 
upper limit on the scattering fraction / sc < 0.25 x 10~ 2 , 
thus Model 6 converged to Model 3. 

In summary, we conclude that Model 3 accounts prop- 
erly for the broad band X-ray continuum measured 
by Chandra. In this model, the nuclear radiation de- 
scribed by a cut-off power-law with a hard spectrum, 
T ~ 1.56, is heavily intrinsically absorbed by a cold gas 
with ATh, z ~ 5 x 10 23 cm~ 2 . The reflected radiation ar- 
rives to the observer unabsorbed, or significantly less ab- 
sorbed, than the intrinsic emission. The upper limit on 
the absorbing column density of the reflected component 
is < 2.5 x 10 22 cixr 2 based on the Model 4 test. 

The reflecting material is neutral and located far from 
the black hole, as suggested by the width and energy of 
the fluorescent iron Ka line, and the ionization parame- 
ter of the reflected component. When this best fit model 
is applied to the Swift/BAT data, we find a normalization 
constant between the two data sets, Cbat = 0.68 ±0.10. 

The right panels in Figure [3] show data-to-model ratios 
(Models 1-3) in the Swift/BAT energy range. The figures 
indicate that the Model 3 scenario is favored and we plot 
the model and its components in Figure [5^,. 

3.2. XMM-Newton data 

In this section we apply the best fitting Chandra model 
(Model 3) to the XMM-Newton data. 

Figure U^i shows the ratio of the XMM-Newton data to 
Model 3 with parameters fixed at the values fitting best 
the Chandra data. A soft excess below 2keV is clearly 
visible, and is due to the diffuse emission of gas surround- 
ing the nucleus collected from a large extraction region 
required by the XMM-Newton PSF. To account for this 
difference we thawed the temperature and normalization 
of the thermal bremsstrahlung component in Model 3 



(kT and Atb, respectively). The fit in the soft band 
improved significantly, although certain residuals indica- 
tive of emission/absorption lines were still present. Since 
we are interested primarily in the continuum AGN emis- 
sion, we do not improve further on the soft (< 2keV) 
residuals and focus on modeling the spectrum at ener- 
gies > 2 keV. The resulting data-to-model ratios are plot- 
ted in Figure |4Jd. The temperature kT is lower than in 
the case of the Chandra spectrum (although consistent 
within the uncertainties) and close to the temperature 
of the thermal component in the Chandra background 
model. Residuals below and above ~ 4keV are still ap- 
parent indicating that the intrinsic spectral properties 
of the nucleus might have changed between the XMM- 
Newton and Chandra observations. The shape of the 
residuals suggests that either the intrinsic absorption col- 
umn varied (model excess below 4 keV) or the normaliza- 
tion of the power-law varied (model deficit above 4keV). 
We test both possibilities. 

First we thaw the intrinsic absorption column (in ad- 
dition to kT and Atb ; Model 3a in Table [2]) and repeat 
the fit. We obtain an increased intrinsic column den- 
sity, ATh, z ~ 8 x 10 23 cm~ 2 , with respect to that in the 
Chandra observation. The intrinsic EW is in agreement 
with that in Chandra Model 3, as expected, while the 
observed EW increases to EW G = 247 ± 10 eV (Tabled 
FigureBh). The fit resulted in CSTAT = 228 (191 d.o.f.). 

Next, we test the hypothesis that the spectral variabil- 
ity is due to the variable normalization of the power-law, 
Nqpl (Model 3b in Table EJ Figure HH), while the intrin- 
sic absorption and reflection parameters remain fixed at 
the values obtained from the analysis of the Chandra 
data. To correctly account for the constant absolute flux 
of the reflected component, we release the condition that 
the normalization of the reflected component is equal to 
the normalization of the power-law. Instead, we fix the 
normalization of the reflected component at the value 
found in the Chandra data. We obtained CSTAT = 209 
(191 d.o.f.). We find that Acpl is lower by a factor 
~ 2 in the XMM-Newton data than in the Chandra 
data. We note, however, that this parameter is subject 
to large errors in the Chandra data. This spectral vari- 
ability scenario results in an increased equivalent width 
of the iron line (observed EW = 280 ± 13 eV, intrinsic 
EWi = 122±6eV) as compared with the Chandra Model 
3. A combined variability scenario (i.e., both Ah, % and 
Acpl allowed to vary; 190 d.o.f.) resulted in parameter 
values and CSTAT consistent with those of Model 3b. 

Finally, we checked the consistency between the XMM- 
Newton and Swift/BAT data. We found C B at = 0.68 ± 
0.10 and C B at = 1.10 ± 0.16 in both cases, Model 3a 
and 3b, respectively. The data-to-model ratios in the 
Swift/BAT energy range are presented in Figures HJ;-d. 

4. DISCUSSION 

We studied the nuclear X-ray emission and X-ray spec- 
tral variability of a low-redshift radio galaxy 4C+29.30 
over a broad 0.5-200 keV X-ray band using our deep 
~ 300ksec Chandra observation, ~ 17ksec of cleaned 
XMM-Newton data, and archival 58-months catalog 
Swift/BAT data. We tested several spectral models, and 
concluded that Model 3 in Table [5] provides the best fit 
to our Chandra data supplemented by the Swift/BAT 
data. In this model the 2-200 keV X-ray emission of the 
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central AGN is decomposed into a strongly intrinsically 
obscured hard cut-off power-law continuum (T ~ 1.56), 
and significantly less absorbed reprocessed emission in- 
cluding a neutral reflection hump and an accompany- 
ing neutral narrow iron Ka line. The soft 0.5-2 keV X- 
ray band is dominated by a thermal diffuse emission, 
and may contain a marginal contribution of a scattered 
power-law component. The same model applied to the 
earlier XMM-Newton data led us to explore two X-ray 
variability scenarios. In the first one, the variability is 
apparent and caused by the variations of the intrinsic col- 
umn density. In the other, the variability arises due to 
changes of the normalization of the power-law emission. 

The hard X-ray spectrum, renewed nuclear radio activ- 
ity (van Breugel et al. 1986; Jamrozy et al. 2007; Liuzzo 
et al. 2009), and weak optical continuum dominated by 
the host galaxy (van Breugel et al. 1986) suggest an 
analogy between the state of the nucleus in 4C+29.30 
and a hard or hard-intermediate spectral state of Galac- 
tic black hole X-ray binaries (GBHs; see e.g., Fender et 
al. 2004; Done et al. 2007). The spectral ages of large 
scale radio structures in 4C+29.30 differ by ~ 100 Myr. 
This timescale, when scaled down by the ratio between 
the black hole mass in 4C+29.30 (~ 10 8 M Q , Section|4~! 
of the present paper; Paper I) and that of a typical GBH 
(~ 10 M Q ), is comparable with the time between consec- 
utive outbursts of GBHs ranging from years to decades 
(for examples of GBHs with at least two documented 
activity periods see e.g., Table 1 in Remillard & Mc- 
Clintock 2006; Figure 2 in Done et al. 2007). Thus, 
it is possible that we observed 4C+29.30 at an early 
stage of its renewed nuclear activity, which could have 
started as recently as several decades ago, judging by 
the spectral age of the youngest nuclear radio structures. 
Interestingly, the age of the youngest radio structure, 
~ 15 yr, found by Liuzzo et al. (2009) is in agreement 
with the suggestion of Jamrozy et al. (2007) that a pow- 
erful outburst took place in 4C+29.30 in the period be- 
tween 1990 and 2005. The X-ray loudness parameter 
(e.g., Zamorani et al. 1981; Avni & Tananbaum 1986), 
a ox ~ 0.3838(^^0 - logF„ jX ), where F u<0 and F v>x 
are rest-frame flux densities at 2500 A and 2keV, respec- 
tively, can be calculated for 4C+29.30 assuming the AGN 
contribution to the optical band as in van Breugel et al. 
(1986), and 2keV flux corrected for absorption (Model 
3). We find a low value, a ox ~ 1, further supporting the 
idea of a hard or hard-intermediate spectral state of the 
source (e.g., Sobolewska et al. 2009; 2011). 

Below we discuss in detail the implications of our mod- 
eling on the nature of the studied AGN and its surround- 
ings. 

4.1. Origin of the X-ray spectral variability 

The three X-ray observations of 4C+29.30 available 
to date, spread over a period of several years, allow us 
to address the question about X-ray variability of the 
source. The early Chandra data (Apr 2001; Gambill 
et al. 2003) were described with a simple power-law 
model (r fixed at 1.69) modified by the Galactic ab- 
sorption (Ah = 4.06 x 10 20 cm -2 , slightly different than 
the value adopted in the present paper), and an intrin- 
sic absorption Ah, z = 4.8^ 4 x 10 23 cm -2 in excellent 
agreement with that found in our new Chandra data, 



7Vh, z = (4.96 ± 0.15) x 10 23 cm -2 . The reported 2- 
lOkeV absorbed flux was 1.1 x 10 -12 ergs -1 cm -2 . Based 
on the Gambill et al. model parameters we evaluated 
that the unabsorbed 2-10 keV flux of the source in April 
2001 was Fc t i — 4.8 x 10 -12 ergs -1 cm -2 , compara- 
ble with that at the time of our Chandra observation, 
F c ,2 = 4.7±5j x 10 -12 ergs -1 cm -2 (Feb 2010). How- 
ever, our deep observation revealed that the spectrum of 
the 4C+29.30 core is more complex due to the soft com- 
ponent and reflection features. The model of Gambill et 
al. (2003) when applied to our Chandra data results in a 
reduction of the 2-10 keV flux, Fc,2, by ~ 6%. It is thus 
possible that Fc, i in Gambill et al. (2003) was under- 
estimated by a few percent. Nevertheless, we conclude 
that the 2-10 keV flux and the properties of the intrin- 
sic absorber were approximately unchanged between the 
two Chandra observations. 

However, a discrepancy between our Chandra and 
XMM-Newton data (Apr 2008) is apparent, as demon- 
strated in Figure |4j We proposed two scenarios for the 
X-ray spectral variability of 4C+29.30. We started with 
the good quality deep Chandra observation to establish 
a base-line spectral model, and we found that the XMM- 
Newton data above 2keV can be explained by varying 
only one spectral parameter: the intrinsic absorbing col- 
umn density, Ah, z , or the normalization of the power- law 
component, Acpl- In the former case (variable Ah, z ), 
the XMM-Newton flux is in agreement with both Chan- 
dra measurements, F^ a = Fc 2 — Fc, ij while the intrin- 
sic column density exceeds that measured with Chandra 
by ~ 60%. In the latter case (variable Acpl) we ob- 
tain F N ,b = (2.5 ±0.1) x 10 -12 ergs -1 cm -2 ~ 0.53F C , 2 . 
Cross-calibration studie£@ show that in general the spec- 
tra measured by XMM-Newton (EPIC/PN) and Chan- 
dra (ACIS-S) agree within 10% in the 0.5-2 keV band, 
and within 20% above 2 keV. Thus, the detected Chan- 
dra and XMM-Newton intrinsic flux difference is too high 
to be fully attributed to the cross-calibration issues, im- 
plying that in this scenario the source may be variable 
on timescales of several years in the 2-10 keV band. 

Variations of the intrinsic column density seem to be a 
common phenomenon in AGN on timescales as short as 
weeks (e.g., Risaliti et al. 2005, 2007, 2009, 2010). On 
the other hand, a non- variable narrow 6.4 keV iron line 
and associated constant neutral reflection not respond- 
ing to the variations of the power-law component, inter- 
preted as the reprocessing of the hard X-rays by a remote 
molecular torus, also seem to be common spectral fea- 
tures in both Type 1 AGN (e.g., MCG-6-30-15, Vaughan 
& Fabian 2004; NGC 4051, Miller et al. 2010, Ponti et 
al. 2006; Mrk 841, Petrucci et al. 2007; Ark 120, Nardini 
et al. 2011; a sample of 14 Seyfert 1 galaxies, Bhayani 
& Nandra 2010) and Type 2 AGN (e.g., a sample of 20 
Seyfert 2 galaxies, Risaliti 2002; MCG-5-23-16, Reeves et 
al. 2006; NGC 2992, Yaqoob et al. 2007, and references 
therein; NGC 4945, Itoh et al. 2008, Marinucci et al. 
2012). 

In our modeling, both variability scenarios resulted in 
satisfactory fits (Table [2]). However, a visual inspection 
of the fit residuals in the XMM-Newton and Swift/BAT 
energy ranges (Figure 2]) seems to favor the scenario 

4 http : //xmm. esac . esa. int /docs/document s/CAL-TN-0052 . ps . gz 
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with variable power-law normalization, -/Vcpl (Model 3b) 
over the scenario with variable intrinsic column density, 
ATh, z (Model 3a). We confirmed that model with vari- 
able -/Vcpl is more likely to be correct than model with 
variable -/Vh> by calculating the respective p- values (Pro- 
tassov et al. 2002) based on the series of 1000 simulated 
data sets (using Sherpa's f ake_pha function) for two test 
cases: (a) variable N-r z (null model) vs. combined vari- 
ability scenario (both -/Vh )Z and -/Vcpl allowed to vary 
(alternative model), and (b) variable Nh jZ (null model) 
vs. combined variability scenario (alternative model). In 
case (a) we found p- value < lx 10~ 3 , which implied that 
model with variable ATjj z is unlikely to be true. Con- 
versely, in case (b) the high p- value = 0.32 suggests that 
the model with variable Nqpl is likely to be true. 

In Figure [SJ: we show a comparison between the best 
fitting Chandra and XMM-Newton total models, and the 
unabsorbed power-law model components in the case of 
the favored scenario with variable -/Vcpl- 

4.2. Iron line and reflection hump 

The first short Chandra observation of 4C+29.30 con- 
tained too few counts to detect the iron Ka line or put 
a limit on its EW (Gambill et al. 2003). Our XMM- 
Newton and Chandra observations are thus the first de- 
tections of the Ka iron line in this object. The line was 
narrow (AE ~ 0.06 keV) and centered at the rest frame 
energy close to 6.4 keV. These line properties indicate 
that it originated in a neutral matter far away from the 
black hole, i.e. outer parts of an accretion disk or a dusty 
torus. 

The intrinsic equivalent width of the iron line, EW;, 
measured in our XMM-Newton and Chandra data is in 
the ~ 70-120 eV range, depending on the adopted vari- 
ability scenario (Table [2]). An iron line originating in 
a lamp-post model due to reprocessing of hard power- 
law X-rays with T — 1.9 by an accretion disk observed 
face-on has an expected EW of ~ 140 eV (e.g., George & 
Fabian 1991). However, this EW is sensitive to a number 
of parameters such as the photon index, viewing angle, 
iron abundance, and optical depth for electron scattering. 
Nandra et al. (2007) give the photon index and inclina- 
tion dependent formulas for the EW based on the fits to 
the Monte Carlo simulations of George & Fabian (1991). 
For r ~ 1.56 resulting from our fits and an assumed in- 
clination of 45°, we obtain the predicted EW of ^170 eV. 
This value should be considered as an upper limit due to 
the unknown disk/X-ray source geometry (the reflection 
amplitude in our data is |Q/27r| ~ 0.3 which may sug- 
gest, e.g., a truncated disk geometry or a patchy corona, 
as opposed to f2/27r = 1 in the lamp-post model). It 
is thus possible that at least some of the photons form- 
ing an iron line with an (intrinsic) EW of 70-120 eV, as 
measured in this paper, originated in an accretion disk. 
This would require that the intrinsic absorber and the 
reflector were two separate media. 

Ghisellini et al. (1994) presented a scenario in which 
the iron line is formed as a result of the reprocessing of 
a power-law continuum by a molecular torus (see also 
Krolik et al. 1994; Levenson et al. 2002). The torus 
scenario for the origin of the iron line provides a com- 
pelling possibility of associating the intrinsic absorber 
with the reflector. The EW of the iron line for tori with 
Nn ~ (5 — 8) x 10 23 cm~ 2 and a half-opening angle of 



30°-45° is predicted in the ~ 80-300 eV range, with the 
lower (higher) values corresponding to the face-on (edge- 
on) view. Ghisellini et al. (1994) assumed an illumi- 
nating power-law with T = 1.9. Levenson et al. (2002) 
estimated that the allowed range of EWs increases by up 
to 5% if the spectrum hardens with r = 1.9 changing to 
T = 1.7 (for intermediate viewing angles i < 65°). The 
spectrum measured in our data is even harder, charac- 
terized with r ~ 1.56, which implies even larger allowed 
EWs. Comparing these predictions with the observed 
iron line EWs in 4C+29.30, we find that both Chandra 
and XMM-Newton measurements, EW„ ~ 170-280 eV, 
lie within the range predicted by the torus origin of the 
iron line. 

In Model 3, we used the pexriv+zgauss model of X- 
ray reprocessing. We are aware of certain issues asso- 
ciated with this approach, which include the iron line 
added to a continuum in an arbitrary way, and the as- 
sumption that the hard X-rays are reflected by an ac- 
cretion disk. We argue instead that the reflector could 
be associated with a molecular torus. Models of hard 
X-ray reprocessing by a toroidal reflector are available 
(Ghisellini et al. 1994; Krolik et al. 1994; Levenson et 
al. 2002; Murphy & Yaqoob 2009). Of these, the tab- 
ulated model mytorus of Murphy & Yaqoob (2009) ac- 
counts self-consistently for all components reprocessed by 
an optically thick torus - the zero-order absorbed com- 
ponent, the scattered component, and associated fluo- 
rescent Fe Ka, Fe Kf3 and Ni Ka emission lines. How- 
ever, the model in its present version considers a torus 
with a fixed half-opening angle of #myt = 60°. This 
model tested on our Chandra data requires almost an 
edge on view, i = 85°, in contradiction with the viewing 
angle constrains based on the jet-to-countcrjct flux ratio. 
The data of 4C+29.30 require that the zero-order com- 
ponent is strongly absorbed, indicating that the torus 
half-opening angle 9 < i < 60°. With viewing angle 
fixed at i — 45°, an additional intrinsic absorber affect- 
ing only the zero-order component of the mytorus model 
is needed in order to obtain a good fit, which breaks the 
consistency of the model. We conclude that the available 
mytorus model does not provide at this time enough flex- 
ibility to account for the case appropriate for our source, 
in which 9 < i < #myt = 60°. 

4.3. 4C+ 29.30 - a 'hidden' AGN 

Recently, a new class of AGN has been identified in the 
Suzaku and XMM-Newton follow-ups of the Swift/BAT 
AGN with N H > 10 23 5 cm- 2 (e.g., Ueda et al. 2007; 
Winter et al. 2008, 2009a,b; Eguchi et al. 2009, 2011; 
Noguchi et al. 2009, 2010). The characteristic feature 
of these AGN is a lower contribution of the scattered 
power- law emission to the soft X-ray band, / sc < 3%, as 
compared to that of typical Seyfert 2 galaxies, / sc = 3- 
10% (e.g., Turner et al. 1997; Guainazzi et al. 2005; 
Winter et al. 2009a; Noguchi et al. 2009, 2010). Winter 
et al. (2009a) found that 24% of their entire Swift/BAT 
selected sample (153 sources) belong to this new class of 
AGN based on the / sc < 3% and low soft-to-hard X-ray 
observed flux ratio, Fo.5_2keV / ^2-iokeV < 0.04. Ueda et 
al. (2007) proposed that the reduced / sc is due to a ge- 
ometrically thick torus covering a significant fraction of 
the solid angle. They predicted that in this new class of 
AGN, the [OKI] lines relative to hard X-ray luminosity 
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would be weaker than in the typical Seyfert 2 galaxies 
because of lower flux arriving to the NLR in the pres- 
ence of a geometrically thick torus. A name 'hidden' or 
'buried' AGN was proposed for this new class, missed in 
previous optically and soft X-ray selected samples. 

Noguchi et al. (2009, 2010) identified hidden AGN 
candidates in the 2XMM Catalogue using X-ray hard- 
ness ratios, HR3 and HR4 (for definition see the 2XMM 
User Guide to the Catalogue). These two hardness ra- 
tios calculated for the XMM-Newton data of 4C+29.30, 
HR3=0.09 and HR4=0.83, comply with the Noguchi et 
al. selection criteria. In our analysis, we found that the 
scattered component was not formally required to model 
the data and that the scattering fraction was consistent 
with zero (Model 6), / sc < 0.25xl0" 2 . This value is quite 
low but in agreement with / sc found previously in several 
other Seyfert 2s (Ueda et al. 2007; Eguchi et al. 2009; 
Noguchi et al. 2009; Winter et al. 2009a, b), and implies 
a torus with half opening angle 8 < 20° (Eguchi et al. 
2009 for / sc < 0.5%). The observed soft-to-hard flux ra- 
tio in 4C+29.30 is F . 5 -2kev/i<2-iukeV ~ 0.005. Thus, 
the XMM-Newton colors, / sc , and soft-to-hard flux ratio 
indicate that 4C+29.30 is a hidden AGN candidate. 

Noguchi et al. (2010) showed that the scattering frac- 
tion / sc correlates with the ratio of the reddening cor- 
rected [OIII] A5007A line luminosity to the 2-10 keV in- 
trinsic luminosity £ — i[oni]/^2-iokcV- Thus, to further 
exploit the intriguing hidden AGN scenario in 4C+29.30, 
we investigated the properties of its [OIII] emission line. 
We used the observed [OIII] flux reported by van Breugel 
et al. (1986). We corrected it for the extinction using 
their observed Hq/H^ luminosity ratio of 5.1 and the re- 
lation of Bassani et al. (1999) with the assumed Galac- 
tic gas to dust ratio, i.e. (H a /H ( g) = 3.1. We obtained 
£ ~ 0.1 in 4C+29.30, which means that our source can be 
placed somewhat above the Noguchi et al. correlation, 
however within the scatter produced by the sources used 
to derive the correlation. 

Base on a sample of six AGN, Eguchi et al. (2009) 
speculated that hidden AGN with / sc < 0.5% have 
high reflection amplitudes S1/27T > 0.8, while those with 
0.5% < / sc < 3% have fi/2vr < 0.8. 4C+20.30 with 
| ft/27r | = 0.32i°;^ and / sc < 0.25% fits in the gap be- 
tween the two groups in the fl/2n vs. / sc diagram (sim- 
ilarly to two other sources in Eguchi et al. 2011), in- 
dicating that the distribution of / sc and fi/27r may be 
continuous rather than split into two distinct groups. 

The soft thermal emission identified in our Chandra 
image contributes majority of the 0.5-2 keV flux. How- 
ever, this component may be easily overlooked in the 
data with lower exposure or spatial resolution. We in- 
vestigated a possible effect of neglecting this component 
in the total model (Model 7 in Tabled]). We found that 
the most striking difference between Model 3 and Model 
7 is the softening of the photon index in the case without 
the thermal component, Tm7 = 2.03 ±0.16. This follows 
from the adjustment of the slope of the scatterred un- 
absorbed component (linked in the model to that of the 
CPL) in order to account for the soft X-ray emission but 
not to overestimate the hard X-rays at the same time. 

We conclude that more detailed studies of the soft 
emission in the hidden AGN complemented by better 
quality hard X-ray data (e.g. NuSTAR) are of particu- 



lar importance to recognizing their X-ray properties and 
geometry of the obscuring torus. 

4.4. Black hole mass 

Noguchi et al. (2010) showed that the scattering frac- 
tion, / sc , anti-correlates with the Eddington luminosity 
ratio, L ho i/L E , where L ho \ ~ 30 x £3-10 keV ( a PPly m g 
a bolometric correction factor of 30 typical for luminous 
AGN; e.g., Vasudevan & Fabian 2007). For 4C+29.30 
the limit we found on the scattering fraction translates 
into Lboi/iE > 0.1, which gives M B h < 1-2 X 10 8 M Q . 

Another black hole mass estimate can be done using 
the empirical relation between Mbh and stellar veloc- 
ity dispersion, cr* (Tremaine et al. 2002). Based on 
the optical spectrum taken from the Sloan Digital Sky 
Survey Data Release 7 (Abazajian et al. 2009) the stel- 
lar velocity dispersion of 4C+29.30 is er* = 207.6 km 
s" 1 (using the STARLIGHT synthesis code of Cid Fer- 
nandes et al. 2005), resulting in the black hole mass 
M BH = 1-6 x 10 8 M Q . 

These two new black hole mass estimates are com- 
patible with the discussion presented in Paper I where 
-Mbh ^ 10 8 M was found based on an argument involv- 
ing energetics of the system. Thus, it is possible that the 
mass of the black hole in 4C+29.30 is close to 10 8 M . 

4.5. Multiwavelength perspective of ^C +29. 30 
4.5.1. Infra-red 

X-ray radiation absorbed and reprocessed by an ob- 
scuring torus would eventually emerge in the infrared 
band. Models of the infrared AGN SEDs under the 
assumption of a clumpy structure of a molecular torus 
(Nenkova et al. 2002; Honig et al. 2006; Nenkova et 
al. 2008a, b) are potentially able to constrain the torus 
geometry in sources observed in IR with a high spatial 
resolution (e.g., Gandhi et al. 2009; Alonso-Herrero et 
al. 2012), with the near-IR 1-5 /jm fluxes being most sen- 
sitive to the inclination effects (Nenkova et al. 2008b). 
4C+29.30 was observed in infrared by IRAS, 2MASS, 
and WISE, of which the two latter surveys sampled the 
desired near-IR band. However, in 4C+29.30 the 2MASS 
measurements are dominated by the host galaxj0, while 
the WISE observations with a spatial resolution 6"-12" 
correspond to our 10" Chandra background region rather 
then the core region. Nevertheless, we use the W2 and 
W3 filters at rest-frame wavelengths Aw2 ~ 4.2 fim and 
Aw3 ~ 11.3 jum, to constrain the viewing angle of our 
source. We calculated respective / = vF^/F^i ratios 
(where Fboi corresponds to Lboi ~ 10 45 ergs _1 ) and 
compared them with the 4.5 /im and 12 /jm models of 
Nenkova et al. (2008b; Figure 10; see also Erratum in 
Nenkova et al. 2010). We obtained f 4.3 Mm ~ 0.06 and 
/n.3/jm — 0.08, which gave i ~ 70° and no constraint 
(our /11.3 pm was too low to intercept the 12 /im model of 
Nenkova et al. 2008b) , in contradiction with the i < 60° 
condition followed from the jet-to-countcrjet flux ratio. 
The /4.3pm and /n.3 pm ratios are a factor 2-6 too low 

5 We estimated the host galaxy contribution to the J, H and K s 
bands using the flux measurement at 5400A being a sum of < 80% 
and > 20% galaxy and AGN emission, respectively (van Breugel 
et al. 1986); the E(B — V) = 0.5 mag reddening possibly due to 
a dust lane (van Breugel et al. 1986); and the lOGyr elliptical 
galaxy template of Buzzoni (2005). 
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to comply with the i < 60° constraint. There are sev- 
eral reasons for this discrepancy to consider, (i) The 
bolometric luminosity of 4C+29.30 might have been un- 
derestimated by a factor of a few due to the uncertainty 
of the applied bolometric correction, (ii) The geometry 
of the obscuring torus may differ from that considered in 
Nenkova et al. (2008b), e.g., the case of a geometrically 
thick torus with half-opening angle < 20° evoked for hid- 
den AGN by Eguchi et al. (2009) would imply a torus 
width parameter a > 70° (see Figure 1 in Nenkova et al. 
2008b for the illustration of a) as opposed to a = 45° as- 
sumed to compute the models, (hi) Finally, the ~ 12 /mi 
WISE flux might be partially affected by the 4C+29.30 
ISM environment (see e.g., Goulding et al. 2012 for ar- 
guments in favor of host galaxy origin of the rest-frame 
9.7 /xm Si-absorption feature in Compton-thick Seyfert 2 
galaxies), e.g., a dust lane covering the central part of 
the nucleus, clearly visible in the archival HST image in 
Figure [H 

We verified the third scenario using the correlation re- 
ported in Gandhi et al. (2009) for a number of Seyfert 2 
galaxies between their 2-10 keV luminosities and 12 fj,m 
luminosities. The IR luminosities in Gandhi et al. (2009) 
follow from sub-arcsecond resolution observations, and 
thus most likely they are representative of the true IR 
nuclear emission. In 4C+29.30 the ~ 12 /im luminosity 
is Li2fj, m ~ 8.3 x 10 43 ergs _1 . This is in agreement with 
the expected Li2 M m ~ 8.8l 1 ; 8 x 10 43 ergs _1 (based on 
our unabsorbed 2-10 keV Chandra flux; Model 3), leav- 
ing scenarios (i) and/or (ii) as plausible causes of the 
discrepancy between the Nenkova et al. (2008b) models 
and 4C+29.30 WISE flux measurements. 

4.5.2. Radio 

4C+29.30 is an intermittent radio galaxy with evi- 
dence for several jet ejection events with the spectral time 
ranging from > 200 Myr (the relic structure), through 
< 100 Myr and < 33 Myr, to ~ 15 yr and ~ 70 yr (Jam- 
rozy et al. 2007; Liuzzo et al. 2009). However, the 
radio activity of 4C+29.30 seems rather mild in compar- 
ison with the powerful radio sources. The total jet power 
of 4C+29.30, Lj ~ 10 42 ergs _1 , constitutes only a small 
fraction of its Eddington luminosity, Lj ~ 10 _4 Le, while 
in the most powerful blazars Lj reaches (O.I-I)Le (e.g., 
Ghisellini et al. 2009). Low-power FR I radio galaxies 
are characterized by much lower values of LjjL-& (e.g., 
~ 10~ 3 in the case of M87), but at the same time they 
accrete at low rates (iboi/^E *C 0.1) probably in a ra- 
diatively inefficient mode. To the contrary, we estimated 
that the Eddington luminosity ratio of our low-power ra- 
dio galaxy is quite substantial, Lboi ^ O.ILe- 

The standard radio loudness parameter is typically de- 
fined as the ratio of 5 GHz radio to B-band (4400A) op- 
tical luminosity spectral densities, R = L u ^/L v ^^ with 
R > 10 characterizing radio-loud AGN (as proposed for 
quasars by Kellermann et al. 1989). In the case of 
4C+29.30', one can calculate 90 < R < 430 based on the 
available archival data (5 GHz total flux given in Liuzzo 
et al. 2009 and 4400A flux extrapolated from the allowed 
range of the 5400A nuclear flux, assuming a opt = 1.2; see 
the discussion in van Breugel et al. 1986). However, this 
value should be taken with caution mainly because of un- 
certainties associated with the attenuation of the optical 



nuclear flux due to the obscured nature of the source, 
and host galaxy contribution known to be particularly 
relevant in the case of low-power Seyfert galaxies (Ho & 
Peng 2001; Sikora et al. 2007). 

We compared the radio properties of 4C+29.30 with 
those of 52 hidden AGN candidates compiled from the 
literature based on Ueda et al. (2007), Eguchi et al. 
(2009, 2011), Noguchi et al. (2009, 2010), and Winter et 
al. (2009a, b). In approximately one-third of them (17 
of 52) the scattered fraction was reported to be < 0.5% 
in at least one publication, that is at a level compara- 
ble with the upper limit estimated in 4C+29.30, imply- 
ing a geometrically thick torus with an opening angle 
<20° (Eguchi et al. 2009). Combination of the available 
archival NVSS (Condon et al. 1998) and NED data al- 
lowed us to calculate the standard radio loudness param- 
eter in the case of 34 sources, including 14 of 17 sources 
with / sc < 0.5%. We found that approximately one-third 
of these sources (11 of 34) could be considered radio- loud 
according to the condition R > 10. Among them we dis- 
tinguished five radio-intermediate sources with R ~ 10- 
20, and six sources with a substantially higher radio loud- 
ness, R ~ (2-27) x 10 3 . We caution, however, that simi- 
larly to the case of 4C+29.30, the calculated values of the 
radio loudness are only approximate due to the obscured 
nature of the hidden AGN. Nonetheless, it seems that 
the hidden nuclei are associated predominantly with the 
radio-quiet AGN. 

Finally, we stress that among the 34 sources with mea- 
sured radio loudness, i?, we found only three (~ 9%) 
with / sc < 0.5% and R > 10 (Mrk 348, 3C 33, and 
NGC 612). Hence, 4C+29.30 appears quite unusual in 
that it shows properties characteristic to both radio- 
loud AGN and hidden AGN with extremely geometri- 
cally thick tori. Based on the present data it cannot be 
determined whether a small fraction of radio-loud AGN 
among the hidden AGN population is due to physical 
constraints or an observational bias. 

5. CONCLUSIONS 

We confirmed that the X-ray emission of the AGN 
in 4C+29.30, a radio source hosted by an elliptical 
galaxy, resembles Seyfert 2 type of activity. We es- 
timated that the mass of the central black hole is of 
the order of 10 8 M Q . We decomposed the broad band 
0.5-200 keV X-ray continuum (Chandra/ X MM- Newton 
supplemented by Swift/JiAT observations) into (i) a 
hard (r = 1.5-1.6) cut-off power-law component atten- 
uated due to a strong intrinsic X-ray absorption with 
column density N-r ~ 5 x 10 23 cm -2 ; (ii) a reflection 
hump (|fi/27r| ~ 0.3) and narrow iron Ka line detected 
for the first time in this source originating due to 
the reprocessing of the AGN power-law emission by 
a distant neutral matter; (iii) a soft thermal emission 
of a diffuse gas clearly detected in our Chandra im- 
age dominating the 0.5-2 keV band; (iv) a marginal 
contribution of a scattered power-law emission to the 
soft X-ray band, / sc < 0.25%. We hypothesized that 
the distant reflector and the obscuring matter can be 
both associated with a molecular torus. Based on / sc , 
the ratio of the observed soft-to-hard X-ray fluxes, and 
the XMM-Newton colors we proposed that 4C+29.30 
belongs to the hidden/buried AGN, a new class that 
is emerging among the Swift/HAT hard X-ray selected 
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AGN. We demonstrated that 4C+29.30, a radio-loud 
source with an intermittent jet activity and / sc below 
0.5%, shows properties characteristic to only ~ 9% of 
the presently identified hidden/buried AGN candidates 
with estimated radio loudness parameter. 
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Fig. 1. — Left: Extraction regions for 4C+29.30 core spectra in the Chandra and XMM-Newton observations (1.25" magenta and 25" 
green circles, respectively), and 10" Chandra background ovcrplotted on top of our deep Chandra image. The ACIS-S Chandra image has 
been filtered and includes only 0.5-7keV photons. The original ACIS pixels were binned to 0.123" and smoothed with Gaussian function 
of o~ — 0.615". Middle: Unbinned Chandra spectrum of the 1.25"core region and folded Model 3 (see Sec. 3.1). Right: XMM-Newton 
spectrum of the 25" source region and folded Model 3b (see Sec. 3.2). 




Fig. 2. — An unsharp mask optical HST image of 4C+29.30 with ovcrplotted Chandra X-ray contours (cyan; Paper I). The crosses 
indicate the optical SDSS identified counterparts. The 1.25" Chandra core extraction region is indicated (black circle). A dust lane across 
the optical nucleus is clearly visible. 
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Fig. 3. — Ratio of the Chandra and Swift/BAT data to different spectral models. Chandra data were binned for the clarity of the plot, 
(a) Model 1 - in the Chandra band clear signatures of iron Ka line are visible around 6kcV (~ 6.4 keV in the rest frame); resulting photon 
index is in disagreement with the Swift/BAT data, (b) Model 2 - the ~ 6 keV residuals disappear after including a Gaussian line in the 
model; the photon index is still too low to agree with the Swift/BAT data, (c) Model 3 with a Gaussian line and reflection fits well both 
the Chandra and Swift/BAT data. 
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Fig. 4. — Ratio of the XMM-Newton and Swift/BAT data to (a) Model 3, fitting best the Chandra data; (b) Model 3 with thawed 
temperature, kT, and normalization, Wxb, of the thermal bremsstrahlung component; (c) Model 3a with thawed -/Vtbi kT, and intrinsic 
absorption column density, Nn iZ ; left: XMM-Newton; right: Swift/BAT; (d) Model 3b, thawed JVtb, kT, and normalization of the CPL, 
iV CPL ; left: XMM-Newton; right: Swift/BAT. 




Fig. 5. — (a) Unfolded spectrum (Chandra and Swift/BAT; crosses) and best fitting model (Model 3; thick solid/black). The model 
components are as follows: hard power-law emission (short-dash/blue) absorbed with an intrinsic column density TVjj z ~ 5 X 10 23 cm - 2 ; 
reflected emission including the ~ 6.4keV (rest frame) iron Ka line (thin solid/magenta); thermal diffuse emission (long-dash/orange); 
Galactic absorption. The Chandra data were grouped for clarity of the plot, (b) Same as (a) for XMM-Newton and Swift/BAT data 
(Model 3b). (c) Solid - best fitting Chandra (black) and XMM-Newton (gray) total models. The enhanced soft diffuse emission in the 
XMM-Newton model is due to larger source extraction region than in the case of Chandra. Short-dash - unabsorbed power-law components 
with variable normalization responsible for the intrinsic X-ray variability of the source. 



TABLE 1 

XSPEC/Sherpa terminology of the spectral models. 



Name 


Descript 


ion 


Model 1° 


phabs 


* 


(zbremss + zphabs*CPL) 


Model 2 


phabs 


* 


(zbremss + zphabs*CPL + zgauss) 


Model 3 


phabs 


* 


(zbremss + zphabs*CPL + zgauss + pexriv) 


Model 4 


phabs 


* 


(zbremss + zphabs*CPL + zgauss + zphabs*pexriv) 


Model 5 b 


phabs 


* 


(zbremss + zpcfabs*CPL + zgauss + zphabs*pexriv) 


Model 6 C 


phabs 


* 


(zbremss + zphabs*CPL + constant*CPL + zgauss + zphabs*pexriv) 


Model 7 C 


phabs 


* 


(zphabs*CPL + constant*CPL + zgauss + zphabs*pexriv) 



NOTE - Models 4-6 converged to Model 3. In Sherpa the individual model names are proceeded by xs. 
a CPL = zhighect * zpowerlw. 
b zpcf abs = / pc *zphabs + (1 — / pc ). 
c constant = / sc . 



TABLE 2 

Results of spectral fits to the Chandra, XMM-Newton and Swift/BAT data. 





Thermal 


Intrinsic 


Cut-off 


power-law 6 


Fc Ka linc c 




Neutral 


Norm." 




CSTAT 




bremsstrahlung 


absorption 


(CPL) 








reflection 










kT 


jV TB 


^H.z 


r 


AfcpL 


E AE 


EW 


EWi 


|Q/2tt| 


Cbat 


Value 


Prob. (d.o.f.) 




keV 




10 23 cm" 2 






keV keV 


eV 


eV 










Chandra 


Bkg 


0.63 ± 0.04 


o nl +U.3U 
°' ul -0.27 




-2.28 ± 0.22 


u.uimo_ 0011 












578 


1 (442) 


Model 1 


6.4+15 


46+ ' 04 


4.08 ±0.29 


0.84+.°;^ 


27+H 










0.20 ±0.03 


1129 


~ 10~ 7 (887) 


Model 2 


5.2+1° 


0.46 ± 0.02 


o qi+0.16 
°- yl -0.09 


0.971°;^ 




6.36 ±0.02 0.07 ±0.03 


192+1* 






0.30 ±0.05 


1045 


~ lO" 4 (884) 


Model 3 


0.7lt°;il 


0.83±g : « 


4.96 ±0.14 


1 56+° 15 
± - OD -0.()4 




6.36 ±0.01 0.06 ±0.03 


173+g 




0.32+HI 


0.68 ±0.10 


1016 


0.002 (883) 


XMM-Newton 


Bkg 


1 62+ 2 ' 24 


i.34^:S 




-0.48 


n+ Ui3 

u - lo -0.08 












1219 


1 (1903) 


Model 3a 


0.60 ±0.07 


8 .5ti:t 


7.88 ±0.26 


1.56* 


103* 


6.36* 0.06* 


247 ± 10 


64 


0.32* 


0.68 ±0.10 


228 


0.03 (191) 


Model 3b 


0.60 ±0.07 




4.96* 


1.56* 


52.0 ±2.7 


6.36* 0.06* 


280 ± 13 


122 ±6 


0.32* 


1.10 ±0.16 


209 


0.30 (191) 



a Normalization of the thermal bremsstrahlung, JVtb, in units 10~ 5 X [3.02 X 10~ 15 /(47rD 2 (1 + z) 2 )] X J n e «idV , where n e is the electron density in cm" 3 , rij is the ion density in 
cm -3 , D is the angular size distance in cm, and z is the rcdshift. 

b Normalization of the cut-off power-law, JVcpLi in units 10 -5 photons keV" 1 cm" 2 s" 1 at 1 keV. 

c The observed (intrinsic) equivalent width EW (EWi) calculated relative to the intrinsically absorbed (unabsorbed) continuum. 
d Normalization constant between the Chandra (XMM-Newton) and Swift/BAT models. 
* Parameter fixed at the Model 3 value. 



